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Abstract. 

An extended study of exclusive vector meson production in ep interactions has been 
performed by the HI and the ZEUS collaborations at the HERA collider. Recent 
f^T) ', measurements are reported and discussed within the framework of the dipole model 

and pQCD. 

O 

OO 

c3 : I INTRODUCTION 

o 

>^ ■ The sharp rise of the electromagnetic proton structure function, F2, toward low 

values of Bjorken x, discovered at HERA [1,2], and the observation of a large 
fraction of diffractive-like events [3,4] are attributed to a large gluon density in the 
proton at very low x values, typically x < 0.01. The rise of F 2 with decreasing x 
can be accommodated by the QCD, DGLAP evolution equations [5] in NLO down 
to momentum transfer squared Q 2 ~ 1 GeV 2 [6,7]. This suggests that perturbative 
effects set in at relatively low values of the interaction scale. However, the validity 
of the DGLAP evolution equation is established through fits to data which involve 
many unknown parameters and therefore it may be doubtful. Exclusive vector 
meson (V) production at high Q 2 has been proposed as an alternative method to 
infer the gluon content of the proton [8]. 

High energy elastic V production in Deep Inelastic Scattering (DIS) may be de- 
scribed, in the rest frame of the proton, by the fallowing sequence of happenings [9] . 
The incoming lepton emits a virtual photon, which subsequently fluctuates into a 
qq pair. The life time of such a quark pair fluctuation is long enough so that it is the 
pair that elastically scatters off the proton and evolves, long after the interaction, 
into a V state. 
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The character of the interaction of the qq pair with the proton depends on the 
transverse momenta of the pair. If the transverse momentum is large, the spatial 
transverse separation between the quarks is small and it forms a color dipole, whose 
interaction with the proton may be calculated perturbatively [10]. The leading 
process is two gluon exchange. If the transverse momentum is small, the color dipole 
is large and perturbative calculations do not apply. In this case the interaction looks 
similar to hadron-hadron elastic scattering and the process should proceed through 
Pomeron exchange as expected from Regge phenomenology [11]. 

The qq wave function of the virtual photon depends on the polarization of the 
virtual photon. For longitudinally polarized photons, small transverse size qq dom- 
inate. The opposite is true for transversely polarized photons. The attractive 
features of elastic V production is that, at high Q 2 , the longitudinal component 
of the virtual photon dominates. The interaction cross section for the latter can 
be, in principle, fully calculated in perturbative QCD. Moreover, for heavy vector 
mesons, like the J/^f or the T, perturbative calculations apply even at Q 2 = 0, as 
the smallness of the dipole is guaranteed by the mass of the quarks. 

Independently of particular calculations [12,13], in the region dominated by per- 
turbative QCD, the following features are predicted: 

• the total 7*p — > Vp cross section, <7 7 * p ^y p , exhibits a steep rise with W, the 
photon-proton center-of-mass energy, that can be parameterized as a ~ W s , 
with 5 increasing with Q 2 ; 

• the Q 2 dependence, which for a longitudinally polarized photon is expected to 
behave like Q~ 6 , becomes milder, more like Q~ 4 , due to the sharp increase of 
the gluon density with Q 2 ; 

• the distribution of t, the four momentum transfer squared at the proton vertex, 
becomes universal, with little or no dependence on W or Q 2 ; 

All these features have been observed at HERA [14-17] in the exclusive produc- 
tion of p°, <p, uj and J/^ mesons. However, the agreement between measurements 
and theory is far from perfect. There are many factors that may spoil the agree- 
ment. First and foremost, the measured cross sections contain the possibly soft 
transverse component of the virtual photon, which has to be modeled. In addition, 
the perturbative calculations have the following uncertainties: 

• the calculation of <x 7 * p ^y p involves the so-called skewed parton distribu- 
tions [18], which are not yet well tested and involve gluon distributions outside 
the range which is constrained by global parton density analyses [19]; 

• higher order corrections have not been fully calculated, therefore the overall 
normalization is uncertain and the scale at which the gluons are probed is not 
known; 

• the fast rise of er 7 * p _>y p implies the presence of the real part of the scattering 
amplitude, which is not fully known; 



• the wave functions of the vector mesons are not fully under control. 



In spite of all these problems, it is generally felt that precise measurements of 
(7 7 * p _ > y p as a function of W, Q 2 , t and the mass of the vector meson, My, with a 
separation into longitudinal and transverse components [20], will help to resolve the 
theoretical uncertainties and ultimately lead to a better understanding of the parton 
distributions in the proton as well as of the dynamics of high energy interactions 
in the presence of a large scale. 

One of the challenges in confronting perturbative QCD calculations with data 
is the ability to establish a region where hard interactions dominate over the soft 
component. The soft component is believed to be well described by Regge phe- 
nomenology, according to which at high energy, the Pomeron exchange dominates 
in the production of diffractive states. The parameters of the Pomeron trajectory 
are known from measurements of total cross sections in hadron-hadron interactions 
and elastic proton-proton measurements. It is usually assumed that the Pomeron 
trajectory is linear in t and has the following form: 

a P (t) = a P (0) + ctpt . (1) 

The parameter 0:^(0) determines the energy behavior of the total cross section, 

atot ~ s^- 1 , (2) 

and a'jp describes the increase of the exponential slope b of the t distribution with 
increasing s. The value of a' p is expected to be inversely proportional to the 
typical transverse momenta squared of the partons participating in the exchanged 
trajectory [21]. Indeed a large value of a' p would suggest the presence of low 
transverse momenta, typical of soft interactions. The latest fits of a P (0) [22] and 
a' p [23] are give: 

a P (0) = 1.096 ±0.003 (3) 
a' F = 0.25 GeV~ 2 . (4) 

These values for the trajectory parameters describe the so called soft Pomeron 
trajectory. The non-universality of a P (0) has been established in DIS, where the 
slope of the rise of the j*p total cross section with W has a pronounced Q 2 depen- 
dence [24,25]. The issue of the a' p dependence on the hardness of the interaction 
may be addressed in the study of exclusive V production at HERA. The value of 
a' p may be determined from the W dependence of b, since b is expected to behave 
as 

W 

b(W)=b + 4a F \n — . (5) 
Wo 

The parameter a' F may also be derived from the W dependence of the differential 
cross section da/dt at fixed t, 



^-{W) = F{t)W^ 2a ^~ 2 \ (6) 

where F(t) is a function of t only. The latter approach has the advantage that no 
assumption needs to be made about the t dependence. The measurements of a pit) 
in exclusive J/ip photoproduction show that, while a P (0) for J/ip is larger the soft 
a P (0) value, the value of a' p for J/ip is smaller [16,26,27] than the soft value. 

Another way to investigate the contribution of hard and soft components in V 
production is by projecting out the interactions induced by the longitudinally and 
transversely polarized virtual photons. Due to the respective wave functions struc- 
ture, small qq configurations are dominantly longitudinal while large configurations 
are dominantly transverse. If one assumes that exclusive V production proceeds 
through s-channel helicity conservation (SCHC), the separation into the longitudi- 
nal and transverse components is possible. At small t, the hypothesis of SCHC has 
been directly tested in the data, and small deviations have been observed [28,29]. 
However, the deviations are small enough not to jeopardize the decomposition of 
the cross sections into the longitudinal and transverse components. 

The angular distributions of the decay products of the vector mesons give access 
to the spin density matrix elements, which are bilinear combinations of the helic- 
ity amplitudes T\ v \ where Ay (A 7 ) is the V (virtual photon) helicity [30]. As t 
increases, a possible change in helicity is more likely due to the change in the V's 
direction and the transfer of the transverse momentum carried by the gluons. The 
following features are expected by pQCD models [13,31,32] in DIS for \t\ ^,Q 2 : 

• a constant ratio of the helicity conserving amplitudes, \T n \ / |T 00 |, with t ; 



a yj \t\ dependence for the ratio of the single helicity flip to the non-flip ampli- 
tudes | Toi| / | Too | and |T\ 1 / |T 00 |; 

a linear t dependence for the ratio of the double flip to the non-flip amplitudes 

|Tl-i| / I Tbo I ; 

the hierarchy 

|T 00 | > |T n | > |T i| > |T 10 | > |T\_ x j . (7) 



The differential cross section, of V production has an exponentially falling 
cross section at low t and a power like behavior at high t. An exponential behavior 
is associated with a Gaussian charge distribution and the slope, b, of the exponent 
is related to the width of the distribution and is a measure of the interaction size. 
In a naive picture for elastic scattering, b has contributions from the size of the qq 
pair which decreases with Q 2 and the constant size for the proton. At high Q 2 , for 
all type of elastic V production, b should reach an asymptotic value equal to the 
proton size. 

The cross section for elastic V photo or electroproduction, with small transfer 
momentum to the proton, has the form [33] 



dt 



F 2 {t)F 2 {t,Q 2 )\A{W^Q 2 



where F p (t) and Fy(t) are the form factors of the proton and the V, respectively, 
that account for the probabilities of their elastic production, and A(W,t) is the 
amplitude of the constituent interaction. The cross section of proton dissociation 
V production does not contain the proton form factor 

= F v (t,Q 2 )\A(W,t,Q 2 )\ 2 . (9) 

Vertex factorization [34] predicts that the ratio of proton dissociation over elastic 
differential cross section is Q 2 independent. The ratio of equation (8) to equation 
(9) is equal to the square of the proton form factor. 



II W DEPENDENCE OF V PRODUCTION 

The dependence of the cross section ay = a 7 (*) p _>y p on the center of mass energy, 
W, is shown in Fig. 1 for photoproduction of various vector mesons [35]. While the 
W dependence of light vector mesons (p to 4>) is oy(W) oc W 022 , as expected from 
Regge phenomenology, the J/ip cross section has a steep rise with W, which is a 
signature of a hard process, as predicted by pQCD. 




1 10 1 o 2 



W„ [GeVJ 

FIGURE 1. The photoproduction cross section as a function of W , for different vector mesons. 
The lines show a W dependence, with S values as indicated. 



A change in the W dependence is also seen when moving from low to high Q 2 , 
as shown in Fig. 2, where the W dependence of a p for various Q 2 values [28] is 



plotted . For each Q 2 value the W dependence of a p is fitted with a form W 5 . The 
parameter 5 is related to the exchanged trajectory a P (t) by 5 = 4(a(< t >) — 1). 
To extract «jp(0) from 5, < \t\ >= 1/6 is taken from measured values [28] and the 
value alp = 0.25 GeV -2 is assigned. The resulting a P (0) is plotted as a function of 
Q 2 in Fig. 3 and a marked increase of 0:^(0) with Q 2 is observed. 
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FIGURE 2. The cross section <7 7 *p_^ pp as 
a function of W, for several values of Q 2 . 
The inner error bars are statistical and the 
full error bars include the systematic errors 
added in quadrature. The lines correspond 
to a fit of the form <r p oc W 5 . 
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FIGURE 3. The Q 2 dependence of the 
intercept a F (0). The inner error bars are 
the statistical and non-correlated system- 
atic uncertainties and the outer error bares 
include the variation of the intercept a F (0) 
when assuming a' F = 0, added in quadra- 
ture. The dashed lines represent the range 
of values obtained for the "soft Pomeron" 
intercept [23]. 



The comparison between the measurements of Oj/^ in photoproduction and var- 
ious pQCD calculations [19,36], based on different gluon distributions, is shown 
in Fig. 4 [37]. The pQCD calculations are known to suffer from large theoretical 
uncertainties, therefore no discrimination between the gluon distributions can be 
made from this comparison. Also shown in Fig. 4 is the result of a calculation [38] 
based on the Golec-Biernat Wiisthoff (GBW) model [39] of the photon, whose pa- 
rameters were fitted to the measurements of the inclusive structure function. The 
GBW model gives a good representation of the data, assuming a double Gaussian 
wave function of the J/ip. 

Although the W dependence of J/ip photoproduction cannot be explained by the 
exchange of a universal soft Pomeron, the data can be fitted to the two-Pomeron 
model [40] as shown in Fig. 5 [41]. The parameters of the established soft and 
the proposed hard Pomeron trajectories are (a(0), a' p ) so ft = (1.08, 0.25 GeV -2 ) 
and (a{ti),a F )hard = (1.418, 0.1 GeV -2 ). The relative contributions of the hard, 
soft and mixing terms of the model are found to vary between 0.1 : 0.5 : 0.4 at 
W = 30 GeV to 0.5 : 0.1 : 0.4 at W = 250 GeV. 

In Fig. 6 the W dependence of ctjm is shown for various Q 2 values [42]. The mea- 
surements are compared with theoretical predictions of Frankfurt et al. (FKS) [43] 
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FIGURE 4. The exclusive J/tp photopro- 
duction cross section as a function of W. 
The inner bars indicate the statistical un- 
certainties, the outer bars are the statisti- 
cal and systematic uncertainties added in 
quadrature. The lines are predictions of 
pQCD calculations. 
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FIGURE 5. The cross sec- 
tion cr("fp — > J/V'p) versus W. The full line 
is the prediction of the two-Pomeron model 
by Donnachie and Landshoff. The separate 
contributions of the hard, soft and mixing 
terms are also indicated. 



using CTEQ4M [44], and of Martin et al. (MRT) [36] using the CTEQ5M [7] 
parametrization. The calculations are consistent with the data within the uncer- 
tainty of the measurements. Also shown in Fig. 7 are the values of S, obtained by 
fitting in each Q 2 , the functional form aj/^ oc W . The results may indicate a 
slight increase in S(Q 2 ) at high Q 2 . 



A Test of SU(4) relation 



The HI collaboration [16] observed that when the cross sections for various Vs, 
weighted by the appropriate S77(4) factors (p° : u : : J /if) = 9 : 1 : 2 : 8) 
were plotted as function of Q 2 + My, at W = 75 GeV, they seem to line up on a 
universal curve, as shown in Fig. 8. However there was a slight indication that the 
J/ij} cross sections lie above the curve (see the insert in Fig. 8). The new precise 
measurements [37,42] of J/ip production do not fit this picture, as shown in Fig. 9, 
where the W dependence of light and heavy Vs, appropriately weighted by the 
S77(4) factor, is plotted in bins of Q 2 + My. The J/ip data lie far above the light 
Vs at all scales. It seems, therefore, that the variable Q 2 + M v is a good variable 
for comparing the light Vs but is not adequate when the J/tp is added to the 
comparison. 
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FIGURE 6. W dependence of the cross 
section a{^*p — > J/ipp) compared with the- 
oretical predictions for Q 2 = 3.1, 6.8 and 
16 GeV 2 . The solid lines indicate the results 
of fits to the function W s . The measure- 
ments from photoproduction are also shown 
for comparison. 
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FIGURE 7. The value of S plotted as a 
function of Q 2 obtained from fits to the form 
cr J/v , oc W 5 . 
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FIGURE 8. The total cross sections for 
elastic V production scaled by SU(5) fac- 
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tors, as a function of (Q 2 
W = 75 GeV. The curve corresponds to a 
fit to the HI and ZEUS p° data. The insert 
shows the deviation of the parameterization 
from the data for to, cf> and J/i/j. 
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FIGURE 9. Comparison of V weighted 
cross section values at fixed Q 2 + M 2 scales, 
as indicated in figure. 



B Measurement of the Pomeron trajectory, a P (t) 

The IP trajectory, a P (t), can be measured by fitting the W dependence of the 
7*p cross section in bins of t, ^\t(W) ~ W 5 , where 5 = A(a P (t) — 1). The 
parameters of the IP trajectory were extracted from measurements of p° [45,46] 
and J/ip production [16,37] and are summarized in table 1. The measurements of 
a P (t) and the fitted trajectories are shown in Fig. 10. and in Fig. 11 

The measured trajectories change with the M v and Q 2 . 



Table 1. Compilation of results obtained for ajp(0) and a' n 



V 


a F (0) 


a' F (GcV- 2 ) 


p°, (Q 2 * o) 


1.096±0.021 


0.125±0.038 


J/1> (HI) 


1.27±0.05 


0.08±0.017 


J/tp (ZEUS) 


1.200±0.009±2:Cio 


0.115±0.018l°;0;i 


p° (ZEUS, DIS) 


1.14±O.Ollg;8i 


0.04±0.07±g;J| 
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FIGURE 10. a F (t), as measured for DIS 
p°, for two Q 2 bins, 2 < Q 2 < 6 and 
6 < Q 2 < 40 GeV 2 , compared to p° pho- 
toproduction results. The solid lines are a 
fit to the data. 
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FIGURE 11. The value of a F (t) as a 
function of t as measured for, J/i/j photo- 
production compared with previous results 
for p° photoproduction. 



Ill ANGULAR DISTRIBUTIONS AND HELICITY 

STUDIES 



The separation of the 7*p cross section into the contribution from longitudinal 
photons, <tl, and transverse photons, <tt, can be performed by measuring the V 



density matrix elements, rf- [30]. These matrix elements are used to parameterize 
the angular distribution W(cos6 h , 4> h) & h ) of the helicity angles. The density matrix 
elements, rfj, are related to the transition amplitudes of longitudinal or transverse 
photons into longitudinal or transverse Vs. The transition amplitudes are denoted 
as Tij, where i = 0(1) is a transverse (longitudinal) photon going into a j = 0(1) 
transverse (longitudinal) V. The relations between the transition amplitudes have 
been predicted by pQCD models [13,31,32]. 

A The ratio of longitudinal to transverse cross sections, 

R = a L /a T 

The most extensive measurements [28,46] of R{Q 2 ) were preformed for p° electro- 
production. The ratio R of cross sections induced by longitudinal polarized virtual 
photons to the transversely polarized ones, (Jl/&t-> is related to the spin density 
matrix elements of p° through 

el-(r§ 4 -A 2 )' 

where is a linear combination of the p° density matrix elements, and A 2 is pro- 
portional to the contribution of the SCH non Conservation amplitude over the total 
amplitude. A small breaking of SCHC has been predicted [13] and measured [28,29] 
(A = 7.9 ± 1.6% [29]). 

The Q 2 and W dependence of R is plotted in Fig. 12 and Fig. 13, respectively. A 
continuous rise of R with Q 2 is observed. The Q 2 dependence of R is well described 
by a functional form R = ±{Q 2 /M 2 ) K with f = 2.16±0.05 and « = 0.74±0.02 [46]. 

On the other hand, there seems to be no W dependence of R for fixed values of Q 2 . 
This indicates that and <jt have the same energy dependence, which is fixed by 
the Q 2 value rather than the longitudinal or transverse configuration. 

B A measurement of the t dependence of the helicity 

structure 

The extracted [47] values of r°Q for p° production are presented in Fig. 14, where 
r 04 K l T oo| 2 +|Toi| 2 K OL j^ OT + ffL ) and N = \ T ..\2_ No signi fi cant var i a tion of 

Tqq with t' — \t — t m in\ is observed. This observation implies that the slopes of the 
exponentially falling t distributions for the transverse and longitudinal s-channel 
helicity conserving amplitudes, T 00 and T n , are very similar. 

The combination rQ + 2rj i 1 is presented in Fig. 15 [47]. Tq is expected [13,31,32] 
to be proportional to the product of the dominant non-flip amplitude T 00 and the 
single flip T 01 amplitude, Tq cx -^-Re(T oT i). The T 01 amplitude is expected to be 
the largest helicity flip amplitude. The matrix element, on the other hand is 
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FIGURE 12. The values of R as a func- 
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FIGURE 14. Measurement of 
function of t'. The curve is a result of fitting 
the data with calculation of [13]. 
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as a function of i'. The curve is a result of 
fitting the data with calculation of [13]. For 
SCHC a null result is expected, independent 
of t'. 



proportional to the non-dominant amplitudes r\ x oc ■^i?e(2TnT 1 ' — 2Ti T 1 L 1 ). The 
strong if dependence of the Tq Q + 2r\ 1 combination is thus attributed mainly to the 
predicted [13,31,32] \Zf dependence of the ratio of Tqi to the non-flip amplitudes. 

The values for r^ Q + 1r\ x are shown in Fig. 16. The matrix element Tq is pro- 
portional to the single flip amplitude T i square, 7q oc — ^ oc (v^) 2 , while the 
matrix element r\ x is expected to have the same t dependence as the double flip 
amplitude Ti_i, rj x oc -^(TuT^^ oc t, therefore the if dependence of the combina- 
tion Tq + 2r\ ± is expected to be linear, up to effects of the single and double-flip 



amplitudes in the denominator N. The combination Tq + 2r;[ 1 is significantly dif- 
ferent from zero and negative, which implies violation of SCHC. The sign of the 
combination gives information on the relative strength of the T iT \ and TuTi_ l 
products of amplitudes, |T 01 | 2 > (TnTjLj), and therefore T 01 > T 1 _ 1 . It confirms 
that the T i amplitude is significantly larger than the double flip amplitude in the 
present kinematic domain. 
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FIGURE 16. Measurement of Tq + 2r;[ 1 as a function of t'. The curve is a result of fitting the 
data with calculation of [13]. For SCHC a null result is expected, independent of t' . 



IV THE DIFFERENTIAL CROSS SECTION AT LOW t 

The differential cross section as function of t has been studied for p° electro- 
production. The values of b(Q 2 ) are extracted from a fit to an exponential form, 
% oc exp(fot), for t < 1 GeV 2 in elastic production (b e i) and for t < 2 GeV 2 in 
proton dissociation p° production (b p d), for different Q 2 regions . The data cover 



a kinematic range of 50 < W < 140 GeV and 2 < Q 2 < 80 GeV 2 (elastic) and 
2 < Q 2 < 50 GeV 2 (proton dissociation). 
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• ep->eYp (Itl < 2 GeV 2 ) 
dcr/dt = A*exp(-bltl) 
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FIGURE 17. The slope b of exclusive and 
proton-dissociative electroproduction of p°. 
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FIGURE 18. Exponential slope b of the 
t dependence for elastic p° production as a 
function of Q 2 . 



The values of b e i and b p d are shown in Fig. 17 as a function of Q 2 . The Q 2 
dependence of b is well described by the following functional forms, 

A/(M 2 (1 + R(Q 2 ))), 
b pd {Q 2 ) + &oo , 

\{Q 2 /M 2 y)\ 

3.46 ±0.14, 
4.21 ±0.12 GeV" 2 . 

One clearly sees the decrease in the value of b e \ with Q 2 , approaching an asymptotic 
value, from fit boo = 4.21 ± 0.12 GeV -2 . The slopes b e \ and b p d are described with 
the same Q 2 dependence, in compliance with the expectations of vertex factoriza- 
tion [34]. It is interesting to note that the ad hoc parametrization of bpd, based 
on the Q 2 dependence of R, describes well the measurements of b p d and b e \. The 
parameterization is motivated by the expected Q 2 dependence of the cross section 
and the relation between the probability for longitudinal configuration and the size 
of the qq pair. 

The measurements of b e i, done by HI [28] (Fig. 18), show the same Q 2 behavior. 
It is also evident that at low Q 2 , the HERA measurements lie systematically above 
the low energy fixed target results. This might indicate shrinkage of the diffractive 
peak, at low Q 2 , as W increases. 

In order to test the factorization hypothesis at the proton vertex [34], the ratio 
of da/dt of the elastic and proton dissociation reactions is studied as function of 
Q 2 , for two fixed values of t, as shown in Fig. 19. The ratio is consistent with being 



M<2 ) 

bel(Q 2 ) 



A 



Q 2 independent. This lends support to the hypothesis that factorization holds at 
the proton vertex for p° electroproduction in the region 2 < Q 2 < 20 GeV 2 and 
|i|<0.4GeV 2 . 
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FIGURE 19. Ratio of proton-dissociative and exclusive p° electroproduction processes as func- 
tion of Q 2 for fixed t values, as given in the figure. The lines are the results of a best fit to a 
constant ratio at each t value. 



V CONCLUSIONS 



It is evident that the hardness of the process in the elastic production of vector 
mesons in ep collisions at high energy, is influenced by both initial (e.g. Q 2 ) and final 
state configurations (e.g. My). This is demonstrated by the change of the energy 
dependence with Q 2 and My, and by the decrease of the b slope with increasing 
Q 2 and My. Longitudinal and transverse photons seem to have the same energy 
dependence at fixed Q 2 , as well as the same t dependence. This indicates that 
the qq pairs initiated by a longitudinal or transverse photon, and contributing 
to the production of a given vector meson, are typically of the same size. This 
size is determined both by the virtualty Q 2 of the fluctuating photon and by the 
probability to form the vector meson final state. 

Although the overall picture is consistent with pQCD predictions, the quantita- 
tive calculations need further improvements. 
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